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By taking into account precursor films accompanying nanodroplets on trapezoidal substrates we show that 
on a mesoscopic level of description one does not observe the phenomenon of liquid-gas-substrate contact 
line pinning at substrate edges. This phenomenon is present in a macroscopic description and leads to non- 
unique contact angles which can take values within a range determined by the so-called Gibbs' criterion. 
Upon increasing the volume of the nanodroplet the apparent contact angle evaluated within the mesoscopic 
approach changes continuously between two limiting values fulfilling Gibbs' criterion while the contact line 
moves smoothly across the edge of the trapezoidal substrate. The spatial extent of the range of positions 
of the contact line, corresponding to the variations of the contact angle between the values given by Gibbs' 
criterion, is of the order of ten fluid particle diameters. 
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I. INTRODUCTION 

Recent progress in device miniaturization has led to 
an increased interest in adsorption of liquids on sub- 
strates structured topographically on the micron- and 
nanoscalei"— . The infiuence of the substrate structure 
on the morphology and on the location of interfaces and 
three-phase contact lines present in such systems is of 
particular interest. Already a hundred years ago Gibbs 
pointed out that an apex-shaped substrate can pin the 
solid-liquid-gas contact line^. If a sessile droplet of fixed 
volume is placed on a planar substrate it forms a spherical 
cap with a contact angle given by the modified Young's 
equation^ 



where tt — is the angle of the apex formed by the sub- 
strate faces (see Fig.[T|). This ambiguity of the local con- 
tact angle at the apex is called Gibbs' condition or Gibbs' 
criterioni^"— . 
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where 7^3, 7^;, and 7 denote the substrate-gas, substrate- 
liquid, and liquid-gas surface tension coefficients, respec- 
tively; r is the line tension coefficient connected with 
the occurrence of the circular three-phase solid-liquid- 
gas contact line of radius R. For macroscopic droplets 
or if the droplet is invariant in one spatial direction and 
forms a ridge, one has i? — 00 and the modified Young's 
equation reduces to the original Young's equation^i'^. For 
a detailed account of the subtleties associated with the 
line tension see Refs.fiol and [Til If the substrate surface 
forms a sharp corner and the three-phase solid-liquid-gas 
contact line is located at its apex the modified Young's 
equation is no longer valid. The corresponding local con- 
tact angle a can take any value within the rangeii^i^ii^rJi 



^ a ^ 



(2) 



cylindricalii2i"— pillars the solid-liquid-gas contact line 
can remain pinned at the corresponding sharp edges of 
the substrates for a range of volumes of the droplets pro- 
vided the contact angle fulfills Gibbs' criterion. This 
fact is widely exploited in the so-called Vapor-Liquid- 
Solid growth process of nanowires made of semiconduc- 
tors such as silicon (Si) or germanium (Ge)^"— . In this 
process a metal sessile droplet is deposited on a substrate 
exposed to the vapor phase of silicon or germanium. The 
semiconductor atoms are absorbed by the metal droplet 
which becomes supersaturated by them. The ensuing ex- 
cess semiconductor material precipitates at the boundary 
of the metal droplet with the substrate, activating the 
growth of a semiconducting nanowire. Quite often gold 
droplets are used as a catalyst. 

Three-phase contact line pinning and Gibbs' criterion 
are also crucial for capillary filling in microchannels pat- 
terned by postsi^ii^i^"— and for dewetting phenomena 
on geometrically corrugated substrates^. In the former 
case, depending on the shape and the height of the posts, 
and on Young's contact angle 9, the liquid front can be 
pinned by the posts so that capillary filling of the mi- 
crochannel might stop. In the dewetting case, the mor- 
phology of the emerging holes is modified by the height 
and the structure of the steps, also due to three-phase 
contact line pinning. 

In addition to the surface and line tension coeffi- 
cients present in Eq. ([T]), the mesoscopic description 
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FIG. 1. Illustration of Gibbs' criterion for a liquid ridge which 
is translationally invariant in the direction normal to the plane 
of the cut shown here. The substrate forms an edge (•) such 
that its surfaces meet at an angle n — tp. Far from the edge 
and in thermal equilibrium the liquid wedge forms Young's 
contact angle 6 with the local substrate surface. If the liquid 
wedge on the right is pushed left it maintains its contact angle 
until its three-phase contact line coincides with the edge of 
the substrate. If it is pushed further, the three-phase contact 
line remains pinned and the liquid wedge increases its angle 
with the horizontal substrate surface to a value a > 9 until 
this angle reaches the value a = 6 + ip, which coincides with 
the case of the liquid wedge being far to the left from the 
edge of the substrate. The angle 6 + with the horizontal 
corresponds to Young's angle relative to the tilted substrate 
surface on the left. Accordingly, if the liquid wedge is pushed 
further to the left it slides down the tilted substrate surface 
keeping its local contact angle 9. The same conclusions are 
reached if the liquid wedge on the far left recedes to the right 
passing the edge of the substrate. The ambiguity of the local 
contact angle oi if the three-phase contact line coincides with 
the edge of substrate vanishes in the limit of a planar surface 
i/j — >■ 0. This macroscopic picture assumes that (in this cut) 
all interfaces are straight, geometric lines. 



of sessile nanodroplets takes into account the effec- 
tive interface potential acting between the substrate- 
Hquid and the liquid-gas interface. In thermal equi- 
Hbrium a nanodroplet with a contact angle less than 
180° is connected with the wetting layer of the liquid 
phase adsorbed at the substrate^"—. The shape and 
stability of such effectively two-dimensional ridges or 
three-dimensional droplets have been discussed in the 
literature^"—. Also the dynamics of nanodroplets was 
examined on substrates structured geometrically by rect- 
angular steps^i^. 

These kinds of mesoscopic studies for planar substrates 
have not yet been extended to the aforementioned apex- 
shaped substrate with an arbitrary angle. For such sys- 
tems, here we focus on how the contact angle changes 
when the three-phase contact line crosses the edge of the 
substrate, and whether on the mesoscale the three-phase 
contact Hnc remains pinned to the edge, as it is the case 
in the macroscopic description (see Fig.[T]). 

In Sec. II we describe the density functional based ef- 
fective interface Hamiltonian which enables us to calcu- 
late the equilibrium shapes of the liquid-gas interface in 
the presence of geometrically structured substrates. In 



Sec. Ill we analyze how the contact angle of the liquid-gas 
interface separating the coexisting liquid and gas phases 
varies upon moving across an apex-shaped substrate. It 
turns out that on the mesoscale the three-phase contact 
line is not pinned to the edge of the substrate and the 
contact angle varies in agreement with Gibbs' criterion. 
The shape of a liquid nanodroplet deposited on a trape- 
zoidal substrate is examined in Sec. IV. Upon increasing 
the volume of this nanodroplet the three-phase contact 
line moves smoothly across the edge of the trapezoidal 
substrate while the apparent contact angle changes con- 
tinuously between two limiting values fulfilling a modified 
Gibbs' criterion. The modification stems from the fact 
that one has to take into account the change of the con- 
tact angle of the nanodroplet with its volume. We show 
that the spatial extent of the region within which the ap- 
parent contact angle changes significantly is of the order 
of ten fluid particle diameters and thus remains meso- 
scopic. We summarize and discuss our results in Sec. V. 

II. MODEL 

In order to determine the effective interface Hamil- 
tonian for an interface separating a liquid-like layer 
adsorbed on a substrate from the bulk gas phase 
we employ classical density functional theory (DFT). 
The corresponding grand canonical density functional 
r2([p(r)]; T, /^) is a function of the temperature T and 
the chemical potential /i, and it is a functional of the the 
spherically symmetric interparticle pair potential ?5)(r) 
and of the external potential Vext (r) encoding the influ- 
ence of the substrate. The interparticle potential w{r) 
is split into a short-ranged repulsive part Whs{r) and an 
attractive part w{r): 

w{r) = Whs{r) + w{r) . (3) 

Two models of the attractive part will be discussed: a 
short-ranged Yukawa-type potential and a long-ranged 
van der Waals potential. 

We adopt a simple random phase approximation for 
the density functional^"—: 

+ 1 y'dVy'dV«;(|r-r'|)p(r)p(r') (4) 

+ Jd^r (Vextir) - fi) p{r) . 

The equilibrium number density profile minimizes 
^l{[p{r)];T, fi). The first term on the rhs represents the 
free energy in the local density approximation of the ref- 
erence fluid interacting via the short-ranged repulsive po- 
tential Wfisif)- The external potential V^xtij) acting on 
a fluid particle located at position r stems from its inter- 
actions with all particles forming the substrate, 

Vextir) = [ dr' psWsi\r-r'\) , (5) 
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where denotes the spatial region occupied by the sub- 
strate with homogeneous number density Ps- As an ap- 
proximation we take p(r) = in that spatial region where 
Vext is repulsive; in the remaining part of space Vext is 
determined by the attractive fluid-substrate interaction 

Ws- 

The thermodynamic state of the fluid is taken to be at 
the bulk liquid- gas coexistence line po{T) and sufficiently 
below the critical point. This implies that the bulk cor- 
relation length is comparable with the diameter a of the 
fluid particle. Under these conditions the nonuniform 
number density profile p(r) can be described within the 
so-called sharp-kink approximation 

PsHkCR, z) = PI e(/(R) -z)+ pg Q{z ~ /(R)) , (6) 

where Q{z) is the Hcavisidc function while pi and pg de- 



note the bulk number densities of the coexisting liquid 
and gas phase, respectively. The local position of the 
liquid-gas interface is described in terms of the Monge 
parametrization z = /(R = {x,y)). Thus by invok- 
ing the sharp-kink approximation we disregard the ac- 
tual smooth variation of the density profile due to ther- 
mal fluctuations and due to the long range of the in- 
teractions governing the system which give rise to so- 
called van der Waals tails^i^. For a finite system the 
density functional V,{[p{r)];T, p) in Eq. evaluated at 
p{r) = pshfe(Rjz), can be systematically decomposed 
into a sum of bulk, surface, line, etc. contributions^"— . 

For a system which is translationally invariant along, 
say, the y-direction, the /-dependent surface contribution 
to the density functional is the sum of two terms: 



(7) 



where Ly is the system size in the invariant direction. The first term in the bracket corresponds to the free energy 
functional per length of a free, fluctuating liquid-gas interface: 



Ax 



dx I dz 



dz' / dyw{x—x',y,z — z'). 



(8) 



The second term describes the effective interaction per length of the liquid-gas interface with the surface s{x) of the 
substrate: 



^znti[s], [./]) =ipi - Pg) 



Ax 



Ax' / dz 



s{x') 



Az' / Ay^piw{x-x',y,z-z') - psWs{x-x',y,z-z')^ 



(9) 



Axuj{x,f{x), [s],Lx) . 



The function uj{x^f{x), [s\,Lx) is the surface density of 
the interaction functional and it is called effective in- 
terface potential. The limit L^: — > oo is taken after 
the appropriate leading terms proportional to Lx are ex- 
tracted from the above expressions (see for example, c.f., 
Eq.dn])). 

For small undulations ^ 1 the expression in the 

bracket in Eq. ^ can be approximated by its local form 
which is called the local effective interface Hamiltonian 
of the system: 

f^Ax[-i^TVfW + uj{xJ{x)M)}- 

(10) 

Within the present approximation the surface tension co- 
efficient of the liquid-gas interface 7 is given by 

1 ~^{Pi- Pgf j Arr^w{r). (11) 



III. LIQUID-GAS INTERFACE CLOSE TO AN 
APEX-SHAPED SUBSTRATE 

In this section we aim at finding the equilibrium shape 
f{x) of the liquid- gas interface which minimizes the effec- 
tive Hamiltonian (Eq. ([TU])) for an apex-shaped substrate 
(see Fig. [2]). The adsorption of a liquid phase at this 
kind of a substrate has been investigated in the context 
of wetting phenomenal^. We shall focus on the ther- 
modynamic states below the wetting temperature at the 
bulk liquid-gas phase coexistence line. We consider con- 
figurations which attain a finite width £0 a-t the far left 
hand side of the apex and detach from the substrate with 
an angle a on the right hand side of the apex (see Fig.[2|). 
We aim at determining the range of accessible angles a 
for such configurations. As a first step, we recall the 
results for a planar substrat o^^'^^i^^" — , for which - for 
specific choices of the fiuid-fluid and the substratc-fiuid 
intcrmolecular pair potentials - one is able to carry out 
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FIG. 2. Schematic profile of a liquid-gas interface at an apex- 
shaped substrate with a characteristic angle (p. The liquid 
layer thickness £o = cos{(fi/2)f{x — >■ — oo) on the left hand 
side of system is finite and the interface detaches from the 
substrate with an angle a on the right hand side of the sys- 
tem. The shapes of the interface and of the substrate surface, 
f{x) and s{x), respectively, are measured relative to the plane 
z = running through the apex horizontally. The system is 
translationally invariant in j/-direction. 



which after one integration leads to 

1 _ ^p(/(a;)) 



1 + f'ixY 
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(16) 



where C is an integration constant. Demanding f{x — > 
— oo) = £o implies f'{x — > — oo) = 0, f"{x — > — oo) = 0, 
and Wp(^o) = 0. According to Eq. ([T5|) the integration 
constant C equals 



C= 1 



^p(^o) 

7 



(17) 



For a profile f{x) diverging linearly for a; — )> oo one 
obtains^ 



cos 9 = 1 



(18) 



the whole analysis analytically. 



Planar substrate 



For a planar substrate the effective interface potential 
bjp{f{x)) = uj{x, f{x), s = [0]) does not depend explicitly 
on X. Accordingly the effective interface Hamiltonian of 
the system (Eq. (|TOl) ) is given by 



dx 



+ ujp{f{x)) - LOp{ap{x))^ , 



(12) 



where on the rhs of Eq. ([T^ the free energy per length 
corresponding to the asymptotic configuration (i.e., for 
|a;| — >■ oo) 



ap{x) = £o + (x — Xd) tan 9 Q{x — Xd) 



(13) 



is subtracted so that M'p[f] is finite for ^ oo. The 
contact angle 9 (unknown a priori) is formed by the 
asymptotes in the limits x — > — oo and x — > oo (see 
Fig.ia): 



9= lim arctan /'(a;) . 



(14) 



The parameter Xd determines the position of the intersec- 
tion of the asymptotes, which is also unknown a priori. 

The equilibrium shape / = f{x) of the liquid-gas in- 
terface minimizing this functional fulfills the equation 



'l + J'ixf 



3/2 



^'piK^)) 



(15) 
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Xd X 
substrate 



FIG. 3. Schematic plot of the equilibrium liquid-gas inter- 
face f{x) which attains a constant value lo for x — oo and 
forms a contact angle 6 with the substrate covered by the 
wetting film of equilibrium thickness io. The parameter Xd 
is the lateral position of the intersection of the asymptotes 
f{x — —oo) and f{x — >■ oo). There is a family of solutions 
f{x — d) with the same free energy which follow from shifting 
the profile f(x) laterally by a constant d. 



For a planar liquid-gas interface corresponding to a 
wetting film on a flat substrate, the substrate-gas sur- 
face tension coefficient equals the equilibrium surface free 
energy density 



Isg = Isi + 7 + a;p(4) 



(19) 



where 7^/ denotes the substrate-liquid surface tension. 
Together with Eq. this renders Young's law 



COS! 



Isg - Isl 
7 



(20) 
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1. Short-ranged forces 

In order to find from Eq. (|15p the explicit expression 
for the equilibrium shape of the liquid-gas interface, we 
choose 



ujp{£) = 47 - ae"'^ + 



'2£ 



(21) 



as a specific model for the effective interface potential, 
where £ is the film thickness divided by the bulk cor- 
relation length ^ in the wetting phase. We take the 
dimensionless amplitude a within the range ^ a 1. 
Within this model the effective interface potential attains 
its minimum at £0 = — ln(a/2) with ajp(^o) ~ —0^7. 
The amplitude a is a unique function of temperature and 
a = corresponds to the transition temperature of con- 
tinuous wetting at which the equilibrium film thickness £0 
diverges and the contact angle 9 vanishes (cos 9 = 1~ 
due to Eq. P^ *). The value a = 1 corresponds to 6* — tt/2. 
In Subsec UlI A H all lengths (e.g., / and x) arc measured 
in units of ^. 

Deriving the above expression (Eq. pip *) for short- 
ranged intermolecular forces requires to go beyond the 
sharp-kink approximation (Eq. ([S])), which corresponds 
to set ting the bulk correlation length equal to zero (see 
Refs. kol [65l and l66l). On the other hand, in the case 
of long-ranged intermolecular forces the sharp-kink ap- 
proximation turns out to be not a severe one. Using 
this approximation in the latter case one obtains the ex- 
act expressions for the coefficients multiplying the two 
leading-order terms in the expansion of the correspond- 
ing effective interface potential in terms of powers of 1 /£ 
(see Ref.lii). 

For weakly varying interfaces, i.e., |/'(a;)| ^ 1. one can 
expand the effective Hamiltonian: 

Jfp[f]^J dx{|/'(a;)2+c.p(/(x))-wp(4)}. (22) 

The equilibrium shape of the interface / = f{x) mini- 
mizes the above functional and satisfies the equation 



which upon integration yields 



!(/'(.)) 



2 _ cpifjx)) 

7 2 



(23) 



(24) 



Here and in the following we omit the ovcrbar indicating 
the equilibrium shape of the interface. The parameter b 
is the first integration constant. (Note that the first term 
on the rhs of Eq. ([M]) can be negative (compare Eq. (|18p ): 
therefore the second term must be positive because the 
Ihs is positive.) A second integration renders the general 
solution of Eq. 



/(x)=ln{^^^[(8a-t-e^(--'^)) -326^]}, (25) 



where d is the second integration constant which shifts 
the position of the liquid-gas interface in the horizontal 
direction (see Fig.|3]). We put on note the property 



f{-x;-b,-d) = f{x;b,d), 



(26) 



which allows one to focus on the case 6^0. The deriva- 
tive f'{x) can be rewritten in the form 



2b ' 



(27) 



According to Eq. ([25]) for finite d and nonzero b the 
interface profile diverges for a; — > 00: 

f{x -^00) = b{x-d) -2 \n{2b) . (28) 
On the other hand 

f{x^-cx,) = -ln^=£o (29) 
(which implies £0 ^ In 2), provided 

b = a\/2. (30) 

Otherwise lim^„j._oo f{x) is either infinite or undeter- 
mined. For b = a\f2 the liquid-gas equilibrium interface 
profile takes the form 



/(x)^4+ln{l + ^^} 



(31) 



and the contact angle 9 fulfills the equation 

tan 9^ f[x^<yo) = aV2 . (32) 

Equation implies cos 6* (1 + 2a^)"^/2 > 1 - 
(compare Eq. (|18p ). Thus for a given a the contact angle 
predicted by Eq. (|22p is smaller than the one predicted by 
the fuh model in Eq. ([12]). Accordingly, for Eq. ([22]) the 
range of angles 9 accessible upon changing the parameter 
< a ^ 1 is [0, 6*0] with 6*0 ~ 55° < 90°. 

The other solution f{x) = - ln(a/2) = £0 of Eq. ^ is 
trivial. It corresponds to & = a^/2 and d = cx) in Eq. ()25p . 



2. Boundary conditions at a finite lateral position of the 
three-phase contact line 

One way to determine the parameters b and d of the 
equilibrium liquid-gas interface profile in Eq. (j25p is to fix 
the value of the function /(xq) = /o and of its derivative 
f'{xo) = /q at a finite position xq. This implies 



and 



-bd 



y'(/^)2-2c.p(/o)/7, 



2bef"{fo + b) -8a 



'bxn 



(33) 



(34) 
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FIG. 4. Dependence of the derivative /g = f'{xo) on the 
height of the interface fo — f{xo) in cases for which the inter- 
face attains a finite value f{x = — oo) = £o for x — oo. Note 
that cJp(/o — ^ oo) = 0, LOp{£o) — — a^7, and ~ — ln(a/2) so 
that /o(/o = £o) =0. The plot corresponds to Eq. ((21) with 
a = 2/e rendering £q = 1. The inset shows the liquid-gas 
configuration for the specific choice a:o = and the height of 
the interface fo — 2. Note that all lengths are measured in 
units of ^ (see Eq. (pTj) ). 



In the previous subsection we checked that only for 
b = a\f2 (Eq. ([30| ) the interface attains a finite value for 
X — >■ — oo. Due to Eq. ([55]) this leads to the relation 



/^ = y^2(a2+c.p(/o)/7)^ 



(35) 



which is depicted in Fig. 21 According to Eq. (PT|) the 
allowed values of /o are bounded from below by £o; i-e-; 
/o ^ 4- 

With h = a\f2 and Eqs. ([Ml) and (l35|) the second pa- 
rameter d is given by the equation 



g^av^d = 4a 



.^»(y^l+^p(/o)/(a^) + l)- 



(36) 



For fixed and /o ^ the shape of the interface which 
attains a finite value for x — >■ — oo is determined uniquely. 
For /o < the function / = f{x) is not defined in 
the whole range x £ (— oo, oo) and thus physically not 
acceptable. 



3. Boundary conditions at the flat asymptote 

If one requires a finite value /(xq) for .tq — > 
implies /o = 4 (Eq. dMI), ./o ^ (Eq. 
(Eq. dSni)), and with Eq. dM]) 



— oo, this 
6 = aV2 



4V2/^e 



-a\/2a;o 



(37) 



The value of the parameter d follows from Eq. ([57]) , but 
it depends on the way in which /g vanishes and xq ap- 
proaches minus infinity. In contrast to the case of fixing 
the height of the liquid-gas interface at a finite lateral po- 
sition xo, in the present case one obtains a whole family 
of solutions, which is given by Eq. (|5T|) and parametrized 
by d. The members of this family of solutions differ only 
by a constant lateral shift. 



B. Apex-shaped substrate 

In this section we investigate the equilibrium shape 
of the liquid-gas interface at an apex-shaped substrate 
(Fig. [2]). The substrate is translationally invariant in the 
direction perpendicular to the plane of the figure and 
the surface of the substrate is described by the function 
s{x) = —\x\ tan((^/2). We assume that the adsorbed hq- 
uid layer attains a finite width £o for a; — >■ — oo and the 
liquid-gas interface detaches from the substrate with an 
angle a on the right hand side of the apex. The value 
of the angle a is not known a priori. If the detachment 
occurs far to the right of the apex, for the geometry in 
Fig. [2] one expects a = 9 — Lp/2. For this shape of the 
substrate even for the short-ranged intermolecular pair 
potentials the effective interface potential cannot be ob- 
tained in an analytical form so that the equilibrium shape 
of the liquid-gas interface has to be determined numeri- 
cally. 

In view of this loss of analytic advantage we now con- 
sider long-ranged interactions as they are realistic for ac- 
tual fluid systems. For the attractive parts of the fluid- 
fluid and substrate-fluid pair potentials we 

takeS,-69 



w(r) 



(0-2 + r2)3 



:(r) 



A. 



((t2 + r2)3 



(38) 



where A> and As > Q are the amplitudes of the inter- 
actions while (7 and as arc related to the molecular sizes 
of the fluid and substrate particles. For this model the 
surface tension coefficient in Eq. (|lip takes the form 



7 



At: 
8^ 



(39) 



In the case of an apex shaped substrate and for the above 
interparticle potentials the effective interface potential 
reads: 



7 



s{x') 



This leads to the disjoining pressure 

duJap{x, z) _ TT 



(^ap{xJix))^{pi-pg) dz dx' dz' dy'^piw{x-x',y',z-z') - psWsix-x',y',z-z')^ . (40) 



d2 



[pi ~ Pa) 



^Pl fr ^\ ^sPs fx Z 



(41) 



where 



i\ap{x 1 z) = QXCX&Vl 



'x^ + z^ + l + x ) tan^ + z 



arctan 



^■\/a;2 + z2 + l+x^ cot ^ + 2: 



cos f 2a;2 sin | - {{x^ + 0^) sin f - z^fM^T^ (1 + cos^ f - a;^ sin^ f ) 



(42) 



[1 + (zcos ^ — xsin ^)2] [1 + (xsin ^ + zcos ^)2] 



Here and in the following we use the fluid-fluid interaction 
parameter a as the unit of length, thus setting tr = 1. 
Upon introducing dimensionless quantities 

2 V piJ Api a 

Eq. dH]) reduces to 
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(44) 



1. Equilibrium shape of the liquid-gas interface 

The equilibrium profile / ~ f{x) minimizes the effec- 
tive Hamiltonian 

j^ap[f] = l dx{j(^^i + r{xr-^i + a',j,{xy 

i^apix, f{x)) - UJapix, aQp(x))| 

(45) 



+ 



of the system where the contributions from the asymp- 
tote (see, c.f.. Fig. [6]) 



^ap ('^) " 
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\x\ tan — -\ , 
' ' 2 cos{(f/2) 



(x — Xd) tana— \xd\ tan 



0{xd - x) 



2 cos{ip/2) 



Q{x-xd) 
(46) 



have been subtracted which renders the integral in 
Eq. (|45p finite. The profile f{x) satisfies the Euler- 
Lagrange equation 



fix) 



[1 + /'(x)2]3/2 



n.p(x,/(x))-4n,p(^,4''^ 



(47) 



Here and in the following we again omit the overbar in- 
dicating the equilibrium shape of the interface. 

Equation (|T7l) is integrated numerically. This is a 
second-order differential equation so that fixing the value 
of the function /o = /(O) and its derivative /g = /'(O) 
at a certain point, say x — 0, leads to a unique so- 
lution, similarly as discussed for the flat substrate (see 
Sec. lHI A2|) . We search for solutions f{x) which attain 
aap{x) for a; — > — cx). There it corresponds to the same 
thickness £q of the liquid layer as for the correspond- 
ing wetting film on a planar substrate (see Fig. [5] and, 
c.f., Fig.|6]). The solutions of Eq. (l47l) . which satisfy this 
boundary condition, will be called g{x). In order to find 
them we proceed as follows: 

1. fix /o and /q at certain values (/o > 0, /q ^ 0); 

2. integrate Eq. (|47p numerically within the range 
X G [ii, L2], where x = Li < and x = L2 > are 
the numerically imposed limits of the system size 
on the left and the right hand side, respectively; 

3. compare f{Li) with aap{Li) and /'(Li) with 
<p{Li) = tan(^/2); 

4. if the differences \f{Li)/aapiLi) — 1| and 
\f'{Li)/a'^p{Li) — 1| are not small enough we 
return to step 1 with a different choice of /q, but 
the same choice of /q. 

Typical values of \Li\ arc in the range of tens of a 
whereas L2 can be very large, e.g., lO^a. The neces- 
sity to consider only relatively small values of \Li\ (as 
compared with L2) is related to the fact that a signifi- 
cantly higher accuracy is needed to solve the differential 
equation in the region where the liquid-gas interface is 
close to the substrate. Due to limited numerical accu- 
racy and due to finite system sizes we are not able to 
find (for a fixed value /o) the value of the derivative /q 
which renders exactly the function g(x) with go = /(O) 
and g'g = /'(O). What can be achieved numerically is to 
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find values /q = /< < g'o and f'^^ = f'^ > g'^ rendering 
solutions which for sufficiently large and negative x fol- 
low the asymptote aap (x) and differ only slightly from it 
in the vicinity of a; = Li, as depicted in Fig. [51 These 
functions are called /<(a;) and f>{x), respectively. The 
values of f{Li) and f{L2) change continuously with /q 
so that the contact angle a = aictan g' (L2) is bounded 
from below by a< — arctan/^(i2) and from above by 
a> = arctan/^(L2), i-C., a< ^ a ^ a>. 



liquid-gas interface detaches from the substrate. This is 
defined as the intersection of the corresponding asymp- 
totes (see Fig. [6]). 




FIG. 5. Schematic plots of the liquid-gas interface at an apex- 
shaped substrate for fixed /o = /(O) and with < g'o (a) 
and /> > g'o (b). In the former case /<(I/i) > aap{Li), and in 
the latter case />(Li) < aap{Li). The dashed line indicates 
the asymptote aap{x) (Eq. (|46}). These two plots correspond 
to case (b) in Fig.jS] 



2. Gibbs' criterion 

In order to access Gibbs' criterion we analyze the de- 
pendence of the results of the procedure described in the 
previous subsection on the choice of the value of /o- As 
in the case of a planar substrate (Sec. lIII A| there ex- 
ists a minimal value of /o, denoted as /q"™, such that 
solutions g{x) of Eq. ^ exist for /o > The 
solution g™™{x) corresponding to /™™ is symmetric: 
5™*" (—a;) = f/™*"(a;). For solutions g{x) correspond- 
ing to /o > /o"™ we define the contact angle a = 
aictang' {L2) and the parameter Xd, where xa fulfills the 
equation s{xd)+(.o/ cos((p/2) = g'{L2){xd - ^2) +5(-^2)- 
The parameter Xd characterizes the position at which the 



z = 




FIG. 6. Schematic shapes of liquid-gas interfaces at an apex- 
shaped substrate. The parameters a and Xd characterize the 
equilibrium liquid-gas interface f{x) and are defined in the 
main text. Various choices for fo = /(O) lead to the cases (a), 
(b), and (c). For — >■ 00 the local contact angle a + ip/2 on 
the far right side approaches the contact angle S on a planar 
substrate. 



Upon increasing fo > /™"' the contact angle a{fo) 
increases and the parameter Xd{fo) decreases, i.e., the 
three-phase contact line approaches the apex. Changing 
the value of fo enables one to plot the dependence of 
the contact angle a on the parameter Xd (Fig. [7]). The 
difference between q;> and «< is so small that the error 
bars of a are not visible on the present scale. 

In the case of a planar substrate the free energies corre- 
sponding to the asymptotic configurations (Eq. (fT3)) ') are 
the same for each interface profile. Therefore the task of 
finding the equilibrium configuration, i.e., the profile with 
the lowest free energy (Eq. P^ ). is posed well. In the 
case of the apex-shaped substrate, for liquid-gas configu- 
rations fulfilling Eq. (|47p with appropriate boundary con- 
ditions, the free energies of the corresponding asymptotic 
configurations, which represent different constraints, are 
different (Eqs. ((45)) and (|46p ). Thus, comparing free ener- 
gies corresponding to different configurations amounts to 
compare free energies characterizing different constraints. 
These free energies as function of Xd can be interpreted 
as the potential of the effective interaction between the 
three-phase contact line and the apex. 

For Xd — > ±00 the local contact angle 
a = aictan g' (L2) tends to its limiting values 9 ^ ip/2 
from below, which are those expected from Gibbs' 
criterion for this geometry (compare Eq. ([2]) which holds 
for the geometry shown in Fig.[T|). For Xd > x*^ the 
local contact angles are, slightly, smaller then 9 — 
(Fig. [7]). The spatial extent of the region within which 
the contact angle a changes significantly can be chosen 
as the region of Xd where 9 — ^ a ^ ^ -I- 0.9 f/2. 

For (y9 = 7r/3 and for parameters of the effective interface 
potential rendering £0 = 2(t and 9 = 7r/4, this width 
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equals d w 11 cr and thus it is mcsoscopic. 
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FIG. 7. Dependence of the contact angle a on the parame- 
ter Xd characterizing the position of the liquid-gas interface 
detachment from the substrate (see Fig.|6]). The horizontal 
dashed lines from bottom to top indicate the angles 6 — ip/2, 

9 - 0.9 (f/2, 9 + 0.9 (p/2, and 6 + ip/2, respectively, where 9 is 
the contact angle on a planar substrate and tt — is the open- 
ing angle of the apex. The quantity d measures the width of 
the region within which the contact angle changes between 
the values a = 6 — 0.9 (p/2 and a = 6 + 0.9 (p/2 and is the 
position above which a < 6 — ip/2. These data correspond to 
Li = -30cr, L2 = lOV, (p = 7r/3, a = 0.5, A = 0.82, and 
p — 0.023 such that the effective interface potential renders 

10 = 2a and 9 = 7r/4. All lengths are measured in units of a. 
The break in slope of a{xd) at a::^; = is caused by the dis- 
continuity in the derivative of the function aap{x) (Eq. (|46|) ) 
which enters the definition of the point Xd. 



IV. SESSILE DROPLETS ON TRAPEZOIDAL 
SUBSTRATES 

A. Planar substrate 

As preparatory work, in this subsection we discuss the 
shapes of interfaces characterizing sessile droplets on pla- 
nar substrates (see Fig.[S]). We assume that the system 
under consideration has a finite extent in x-direction, 
^ = L, and is translationally invariant in the 
y-direction. The shape of the interface is described by a 
function / = f{x). The total volume Vtot/Ly of liquid in 
the system per length Ly is fixed. Ly denotes the size of 
the system in y direction. 

Within our mcsoscopic description for weakly varying 
liquid-gas interfaces the corresponding effective Hamilto- 
nian is given by (compare Eq. ([H])) 

/L 
dx[lf'ixr+ujpif{x))}, (48) 

where the limits of the x-integration reflect the finite ex- 
tent of the system. The equilibrium shape f{x) of the 



z = 





substrate 




-* *- 



2L 



FIG. 8. Schematic equilibrium shape of a ridgelike liquid 
nanodroplet deposited on a planar substrate. The system is 
translationally invariant in y-direction and has a finite lateral 
extent \x\ ^ Lx = L. 



interface minimizes the functional 

^;[f]=^^p[f]-XVtot/Ly, (49) 
where A is a Lagrange multiplier, and 

Vtot^Ly j Axf{x). (50) 

The equilibrium profile satisfies the differential equation 

7/"(x) =..;(/») -A. (51) 

In the following we again omit the overbar denoting the 
equilibrium configuration. In x-direction we impose Neu- 
mann and periodic boundary conditions: 

n-L) = f'{L) = Q, f{^L) = f{L) = fL, (52) 

with the thickness not fixed a priori. The conditions 
f'{—L) = f'{L) = can be realized by vertical sidewalk 
at a; = ±L exhibiting a contact angle of 90°. 
Integrating Eq. (|5ip renders 

y'{xf^uj^{f{x))^\f{x)~C, (53) 

where the integration constant C is determined by the 
boundary conditions: 

C^OJp{fL)-\fL, (54) 

which leads to 

y'{xf = ujpifix)) - upih) - A(/(x) - h) . (55) 

We examine effective interface potentials ujp{£) with a 
minimum sit £ = £q, uj{£ — ^ 00) = 0~, and one inflection 
point at £ = £1 > 4 (see Fig.|n|). 

We search for solutions /(x) ^ /l, such that 
f'(x — ±L) > and that there is one Xq G {~L,L) for 
which f'{xo) = 0; we denote the maximum value of the 
function f{x) as /o = f{xo). As a result one obtains from 
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Eq. (|51|) the relation A ^ Ld'p{fL) and from Eq. ([55]) one 

has 



A 



/o ^ /l 



From the structure of Wp(/) one infers A > 0. Thus 
the slope of the line connecting the points (/l,Wp(/l)) 
and ifo,oJp{fo)) is equal to the Lagrange multiplier A. In 
order to fulfill the condition A ^ oj'p{fL) this line segment 
must be located below the graph of the effective interface 
potential (see Fig. [9]). This implies two restrictions: one 
for the thickness of the liquid layer at the boundaries 
/l > ^Oj the other for the maximum value of the 
function /o > ii- 



UJp{i) 




FIG. 9. Schematic plot of the effective interface potential 
considered here in the context of discussing the shape of small 
droplets (Eqs. (|55|l and (|57|) ). The slope of the dotted line 
connecting the points {fo,^p{fo)) and (fL,i^p{fL)) equals the 
Lagrange multiplier A (see Eq. (|56p '). 



Integrating Eq. (|55p one obtains 



fo 

dz 



ujp{z) - ujpijL) ~ - h) 



-1/2 



(57) 



On the other hand, by rearranging the limits of inte- 
gration on the rhs of Eq. ((57)) it can be expressed as 
L ^ xq which renders xq = 0. In addition, the func- 
tion f{~x) also fulfills Eq. ([55]) so that in the following 
we consider functions f{x) which are symmetric with re- 
spect to a; = 0. The excess volume 

Ve. = Vtot - 2L Ly fL (58) 

of the adsorbed liquid is given by (see Eq. ([55]) ) 



dx 



fix) - h 



^-V2Ly^ / d 



fo 



fL Vwp(z) - Wp(/l) - Xiz - /l) 

(59) 



By combining Eqs. - ([5^]) with a given volume 
Vtot and a size L one determines the quantities /l, /q, 



A; integrating Eq. (|55|) gives the shape of the equilibrium 
liquid-gas interface in the form (for a; > 0) 



(56) x{f) = L 



dz 



1 



fL V^pi^) - ^pUl) - A(z - /l) 



(60) 



In the following, for the surface tension (Eq. (HI])) and 
for the effective interface potential we adopt the expres- 
sions following from the sharp-kink approximation for 
the density functional using the fluid-fiuid and substrate- 
fluid pair potentials given by Eq.()38|): 



"^pl*) = - ^P\-] ~ ^^P - — 



(61) 



where Cjp{t) = ^arctan(l/£) — 1, and the three parame- 
ters /5, tTs , A are given by Eq. (|43)) . Instead of using the 
parameters A and p, we use the width £o minimizing the 
effective interface potential u!p{t) and the Young contact 
angle given by cos 9 = 1 + a;p(£o)/7. For a given Gg the 
relation between {A,p) and {(-0,9) is unique. 



1. Stability of the droplet 

In addition to the droplet-like solution, Eq. ([55]) has the 
trivial, flat solution f{x) = with zero excess volume 
and nonzero total volume. We flx the lateral size L and 
check which of these solutions has the lower free energy 
(Eq. (|48pL and thus corresponds to the stable interface 
configuration. 

For a finite system with a prescribed fixed L the droplet 
solution does not exist for arbitrary fi^ > Iq. This is 
caused by the boundary conditions f'{—L) = f'{L) ~ 0. 
If /i is sufficiently small, i.e., /l < flW (where the 
threshold value /£ depends also on the parameters of 
the effective interface potential), there is no /o which, 
according to Eqs. ([ST)) and ([56)) . would correspond to the 
prescribed fixed L fFig.irU]). For > ,fl{L) there are 
two values of fo which correspond to the same fixed L 
fFig.irU]). As expected intuitively, one can show that the 
configuration corresponding to the smaller value of /o has 
always the lower free energy. 

On the other hand, for a given value of L, the quantity 
/l is bounded also from above: /£ ^ /l =^ ^i, rendering 
upper (F™t''^(L)) and lower (\/™"(i)) bounds for the 
total volume for which the equilibrium droplet configu- 
ration exists (Fig. fTTj) . One can check that for > £i 
the droplet configuration has a larger free energy than 
the flat fllm configuration with the same total volume. 
The corresponding phase diagram is presented in Fig.ITT] 

It is worth noticing that there exists a smallest value 
Lmin below which the droplet configuration cannot ex- 
ist. Upon increasing Vtot for L > Lmin the droplet con- 
figuration forms discontinuously at Vtot = Vtot^'i^) and 
ceases to exist for Vtot = ^^t^°^(-^^): also in a discontin- 
uous way. The transition values V^^"^ and VJ™" corre- 
spond to /i being equal to /£ and to £i , respectively. For 
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FIG. 10. Dependence of L on /o, according to Eqs. (|57p 
and (|56p. for three different choices of /z,: (a) /l = 2.1, 
(b) /i = /* = 2.119, (c) /l = 2.15. For each choice of fL, 
fo is larger than a minimal value /o (/l) at which L diverges 
as L{fo — >■ /o) ~ — log(/o — fo). The surface tension co- 
efficient and the potential parameters are chosen such that 
as = 0.5, lo = 2, and 9 = ir/A. If one prescribes the width 
2L of the system (e.g., 2L = 60, dashed horizontal line), de- 
pending on the choice of /z, there are two (c), one (b), or no 
corresponding values for fo. For i = 30 the threshold value 
is /£ — 2.119 (b). All lengths are given in units of a. 
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FIG. 11. Morphology phase diagram in terms of the variables 
(L, Vtot) displaying the coexistence lines between the phases 
with a droplet and a flat film configuration, respectively. The 
upper {Vtot = Vt'SriL)) and the lower {Vtot = Vi^r{L)) solid 
lines correspond to the droplet configurations with /i = /£ 
and /l = ^1, respectively. The dashed straight line shows 
the volume of the flat film configuration with the thickness 
/i = l\. Lmin denotes the smallest lateral width of the sys- 
tem for which the droplet configuration can exist. The surface 
tension coefficient and the potential parameters are chosen 
such, that (Ts = 0.5, £o ~ 2, and 6 = 7r/4. Note that all 
lengths are measured in units of a and Vtot is measured in 



L < Lmin the system extension 2L turns out to be too 
small to accommodate the droplet and to simultaneously 
fulfill the boundary conditions /'(— L) = f'{L) = 0. For 
i — > oo the minimal volume increases linearly as function 
= 2L€i (so that for large L in Fig.Othe lower 
bound of the droplet phase approaches the red dashed 
line) and the first-order character of the transition be- 
tween the flat film and the droplet configuration weakens 
and becomes continuous at L = oo. 



2. Macroscopic system 

In this section we recall the analyses of a liquid droplet 
adsorbed at a fiat, unbounded substrate which can spread 
over the whole macroscopic extension of the substrate 
{L — >■ oo)^"—'^. If one places a droplet onto an in- 
finitely extended, fiat liquid film, after some time the 
droplet will disappear into the film without changing the 
thickness of the latter, because the droplet volume is van- 
ishingly small compared with the total liquid volume. 
Accordingly, in theory here we fix the excess volume of 
the droplet above the flat film. In practice, providing 
an experimental setup which on one hand mimics infi- 
nite substrate extensions and on the other hand allows 
for the persistence of a droplet configuration seems to be 
rather challenging. 

The effective Hamiltonian for such a liquid-gas inter- 



face / ~ f{x) with a rfrop is given by 

J-^ ^ ^ ' (62) 

+ u;p(/(a:)) - Wp(£oo)| , 

where £oo is the, a priori unknown, height of the equi- 
librium liquid-gas interface at infinity. The equilibrium 
profile / = /(x) minimizes the functional 



(63) 



where A is a Lagrange multiplier, and the excess volume 
is given by 



dx 



fix) -to 



(64) 



For X ±oo one has f'{x) = f"{x) = so that there the 
interface profile attains a certain finite height too ^ ^o- 
We recall that minimizes the effective interface poten- 
tial, i.e., oj'pito) = 0, and it equals the equilibrium thick- 
ness of a planar liquid film adsorbed at a flat substrate 
in a grand canonical system in contact with a reservoir. 

The equilibrium profile f{x) (here and in the following 
we omit the overbar indicating the equilibrium profile) 
satisflcs the equation (compare Eq. (|5ip ) 



f"{x) 



(l + /'(z)2)3/2 



c^;(/(^))-A, 



(65) 
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which gives 
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-LUp{fix))+Xfix)+C. (66) 



The integration constant C and the Lagrange multipher 
A follow from the boundary conditions at infinity: 



and 



(67) 



(68) 



The equation for the liquid-gas interface configuration 
reads: 

VI + f (x)^ (69) 

+ Oj'pi£oo)if{x)-eoo). 

We search for solutions f{x) ^ £oo- The values of the 
function f{x) describing the equilibrium interface config- 
uration lie in the range [£oo, fo] (Fig. [T2|) . The boundaries 
£oo and /g of this interval are the solutions of equation 
iOpi£) = io'pi£oo){i~ioo)+i^pii^)- The value fo = /(xq) 
is the maximum value of the function f{x). It is a func- 
tion of £ao given by (Eq. ([69l) ) 



^;(^oo) 



(70) 



This means that the line tangent to the curve ujp = ujp{£) 
at the point £ = £oo intersects the curve tOp = u!p{£) 
again at the point £ = fo (Fig. [T^)> ^. The function 
f{x) is symmetric with respect to the x = xo, i.e., 
f{xo + x) = f{xo - x). 




FIG. 12. Schematic plot of the effective interface po- 
tential. The values of the function f{x) describing 
the equilibrium interface configuration lie in the range 
f{x) € [ioojfo]', £oo and /o are the solutions of the equation 
LOp{£) = LOp{£oo){£ -£oc,)+ tOp{£oo) (dotted line). The widths 
£o and £i denote the abscissa of the minimum and of the in- 
flection point of the effective interface potential, respectively. 



inflection point, ujp{£i) = 0, the thickness £oo of the liquid 
layer at infinity is restricted to £o < £oo ^ (-i and thus 
/o ^ ^1 . (These inequalities follow from the fact that the 
dotted line in Fig. [T^] is tangent to ujp{£) dX £ = £oo.) For 
(■oo = ii one has the flat solution f{x) = £oo- 

For weakly varying liquid-gas interfaces (|/'(a;)| ^ 1) 
Eq. (j69|) reduces to (compare Eq. (j55p in conjunction with 

Eq. dnzi)) 



^f{xf =u:p{f{x))~u:p{£^) 



^'p{too){f{x)~£oo). 



(71) 



The excess volume of the adsorbed liquid can be 
expressed as (compare Eq. (|59p in conjunction with 
Eq.dnil)) 



dx 



f{x)-£o 



Lyv/27 / dz 



fo 



z-£. 



^Ju}p{z)-U}p{£oo)~Uj'p{£oo)iz- £oo) 

(72) 

Combining Eqs. (1701) and (|72p with a given expression 
LOp{£) for the effective interface potential and a given ex- 
cess volume Vex one is able to determine the quanti- 
ties £aD, fo'i integrating Eq. ([7T|) gives the shape of the 
equilibrium liquid-gas interface configuration as (com- 
pare Eq. (|60| in conjunction with Eq. (|67l) ') 



xif)^xo- 



fo 

dz 



1 



/ ^UJp{z)-UJp{£oo)-i0'p{£oo){z~£oo) 



(73) 



for x > Xq. Due to the translational invariance of the 
substrate the position xo is finite and arbitrary. The 
solution of the second order differential equation for the 
equilibrium shape of the interface (Eq. (|65|)) contains two 
integration constants: xo and £oo- The latter one can 
be determined from the fixed excess volume constraint 
(Eq. (|72p ) once fo is expressed in terms of ioo by using 

The typical shape of the equilibrium liquid-gas inter- 
face is shown in Fig.[T31 The function f{x) is symmetric 
around the position xo of the maximum. The excess vol- 
ume Vex of the drop uniquely determines £oo and fo (see 
Eqs. and ([70])). For Vex they both reach, with 
vanishing slope, the position £i of the inflection point 
of the effective interface potential (sec Fig. (fT4|)). For 
Vex ^ oo the maximal height fo grows without limit and 
approaches io. 



3. Contact angles 



For an effective interface potential ujp{£), as shown in 
Fig. [121 with ujp{£ 0) > 0, u}p{£ oo) = 0~, and one 



For nanodroplcts the definition of the contact angle 
requires more care than for macroscopic drops. We define 
the contact angle of the nanodroplct as follows fFig. [T51) : 
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FIG. 13. Shape f{x) of the ridgehke equihbrium hquid nan- 
odroplet and the fitted arc of a circle ar{x). The contact 
angle for the droplet is denoted by 9d. The effective interface 
potential parameters are as = 0.5, 9 — it/ 4:, and £o — 2; xq is 
chosen to be 0. The excess volume is such that £oo = 2.1 and 
fo = 8.61. All lengths are measured in units of a. 



3. define the contact angle of the droplet as 



9d = arctana^(a;£i) = arctan 



Xd 



(74) 



One can show, that the cosine of the contact angle of 
the nanodroplet as defined above is given by 



cos 9d 



1 



7 



1 



^(/o) 



(75) 



With decreasing thickness £oo at infinity (i.e., increasing 
excess volume of the droplet, Fig.[T4| the height fo at the 
center and also the contact angle increases^. According 
to Fig.[T31 for Vex — >■ oo one has £oc — >■ and fo — >■ oo 
so that u!p{fo) and uj'p{fo) vanish. Thus, as expected, for 
large droplets the contact angle 9d reaches Young's an- 
gle 9 (see Eq. (|18p and Fig. [T51) . This size dependence of 
9d must be taken into account while investigating Gibbs' 
criterion for a sessile nanodroplet deposited on a trape- 
zoidal substrate (see below). 




FIG. 14. Dependence of the thickness £ao (lower curve) and 
of the maximal height fo (upper curve) on the excess volume 
Vex for the parameters as — 0.5, 9 = 7r/4, and £o — 2. The 
values of £o and £i — 2.65 (compare Fie:. ll2p are marked by 
dashed horizontal lines. All lengths are measured in units of 
a and Vex in units of L^a^ . 




FIG. 15. Dependence of the contact angle 9d of the 
liquid droplet on the thickness at infinity £az for the 
same potential parameters as the ones used in Fig. 1141 in 
particular £q = 2 and l\ = 2.65. The macroscopic limit 
9 = tt/A: is marked by the horizontal dashed line, and 



The limiting value Ud 

N3/2 



e<i = 



vd{l^^lo)-9^^{l^-£^ 

is approached as 9d{£oa — >■ £i) ~ {£i — £ozY^ ^ ■ lengths are 
measured in units of a. There is a one-to-one correspondence 
between £oo and Vex (see Fig. ll4p . 



B. Trapezoidal substrate 



1. find the arc of a circle ar(x) = \/ Iff ~ — i? + /o 
with the same curvature l/i? as the curvature of 
the liquid-gas interface at its maximal height; 

2. find the intersection point Xd di z ^ a^{x) and 



1. Macroscopic description 

In this section we analyze the shape of a droplet of 
a fixed volume which is deposited on a trapezoidal sub- 
strate characterized by the angle (p and by the width 26 
of the planar basis fFig. fTB]) . The system is taken to be 
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translationally invariant in the j/-dircction. The surface 
free energy ^ of the droplet has the foUowing form: 



(76) 



where A, A^i^ and A^g denote the areas of the liquid- 
gas, solid-liquid, and solid-gas interfaces with the cor- 
responding surface tension coefficients 7, 7s;, and 7sg. 
Within this macroscopic level of description the effective 
interaction between the substrate-liquid and the liquid- 
gas interfaces is not taken into account. We focus on 
the case that the droplet is deposited symmetrically on 
the substrate. Moreover, we restrict our analysis to sit- 
uations in which Young's local contact angles Q are re- 
stricted to 6* < 1^/2 — ip. so that even over the sides of 
the trapezoid the liquid-gas interface can be described 
by a single- valued function / = f{x), where the x-axis 
is parallel to the planar basis of the trapezoid. In order 
to describe liquid-gas interfaces with overhangs another 
parametrization is needed, e.g. by the arc length of the 
interface. But then the density functional (Eq. ([7])) has 
a much less transparent form. 

The equilibrium shape of the liquid-gas interface, 
which minimizes the free energy in Eq. (|76p , forms the 
cap of a cylindrical ridge. Depending on the volume V of 
the liquid drop one of three distinct types of configura- 
tions occurs (Fig.[T6|): (/) the area of the substrate-liquid 
interface is smaller than the area of the horizontal basis of 
the substrate and the apparent contact angle is equal to 
Young's angle formed with a horizontal, planar surface; 
(//) the area of the substrate-liquid interface coincides 
with that of the horizontal basis with the three-phase 
contact line pinned to the edge and with the apparent 
contact angle a formed with the horizontal basis in the 
range 6 ^ a + ip; {III) the area of the substrate- 
liquid interface exceeds the one of the horizontal basis 
and the apparent contact angle formed with the tilted 
side of the trapezoid is again Young's angle. 




FIG. 16. Three possible types of cylindrical cap-like configu- 
rations I, //, and /// of a sessile liquid ridge on a trapezoidal 
substrate characterized by the angle ip and the size 2b of the 
planar basis. The system is translationally invariant in the 
direction normal to the figure shown. 



For configurations / and /// the constant radius R of 
curvature of the equilibrium liquid-gas interface and its 
maximal height h above the horizontal basis are given by 



R]{V) 



and 
and by 

and 
hniiV) 



hi[V)^Ri{V){l-cose), 



V/Ly — b'^ tan ip 



(78) 



9 + (fi — sm{6 + ip) cos{0 + (fi) — sin" (6 + (fi) tamp 

(79) 



- 1 ^ cos(6' + (p) — sm{9 + (p) tan Lp 

+ b tan ip . 



(80) 

Ly is the spatial extension of the system in the invariant 
y-direction. 

In configuration // the apparent contact angle a is not 
fixed by materials properties. It is not given by Young's 
equation, but depends on the volume of the sessile droplet 
and is determined implicitly by the equation 



V/iLyb^) 



sin^ a 



cot a . 



(81) 



Equation ((8T|) states that the corresponding section of a 
circle has the area V/Ly. In this case the constant radius 
R of curvature of the interface and its maximal height 
are given by 



Rii{V) 



sma 



and 



hiiiV)^RiiiV)il~cosa), 



(82) 



(83) 



— sni f cos ( 



(77) 



respectively. The dependence of the radius R and of 
the maximal height h of the interface on the volume are 
shown in Fig. 1171 The volumes Vi and V2 are the limit- 
ing values for configuration //, and can be calculated by 
replacing in Eq. ([8T|) a by ^ and 9 + ip, respectively. We 
emphasize that the radius R of the interface is a decreas- 
ing function of the volume for configuration // but an 
increasing function otherwise, regardless of the angles 9 
and If. On the other hand the height h of the droplet is 
an increasing function of the volume in configurations / 
and //. 

For configuration /// the height increases with vol- 
ume ioY > Lp (which is compatible with the constraint 
9 < tt/2 — ip for ip < tt/4) and it decreases with volume 
for 9 < (p. For 9 = tp the height remains constant, i.e., 
it is volume independent. For fixed angle ip there is a 
minimal value 0min(</') < below which there is no one- 
drop solution in configuration ///. The angle 0min(<rf5) is 
the zero of the denominator on the rhs of Eq. ([75]) . For 
^'min(</5) ^ 9 ^ ip the volume of the droplet in configu- 
ration /// is bounded from above by a maximal volume 
¥^■{9, (fi), so that for such a volume the liquid-gas inter- 
face touches the edges of the substrate and the droplet 
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splits into three parts. The radius i?™ of curvature cor- 
responding to V = V™{9, If) is given by 

The analysis of morphological phase transitions between 
distinct sessile droplet configurations on a trapezoidal 
substrate, at which one droplet splits into two or three 
droplets, is left for future research. This has been al- 
ready investigated for droplets deposited on axisymmet- 
ric pillar-like substrates2ii2^. 




= V/{Lyb^) 



III 




with and the effective interface potential between the 
substrate-liquid and the liquid-gas interfaces. The dis- 
joining pressure for the trapezoidal substrate (Fig. [TO]) 
is the difference of the disjoining pressures Ilap{x, z, if) 
(Eq. ([4T|l ) corresponding to two apex-shaped substrates 
(Fig.® 

51 =|(a;,y,z) € : z < {x + b)tainp A z < o| , (85) 

52 =^ix, y, z) e : z < A z > ~{x - b) tamp^ , 

(86) 

for which the characteristic angles arc given hy (pi ~ (p 
and ip2=TT — if, respectively, with x = as the center of 
the trapezoidal basis (see Sec. lIII[) . Thus the disjoining 
pressure stemming from the trapezoidal substrate is 

lltrapix,Z;b,ip) = Ilap{xi,Zi,ipi) - flap (^2 , ^2 , ^2 ) 

(87) 



where 



and 



cos ^ 
-sin^ 



sin^^ 
cos^^ 



cos ^ sin ^ 



fx + b 



- b 



(89) 



are the coordinates corresponding to the above apex- 
shaped substrates 5*1 and S2, respectively. 

For the attractive parts of the fluid- fluid and substrate- 
fluid pair potentials of the van der Waals type (Eq. ([38])). 
the disjoining pressure is positive near the substrate, has 
two saddle points, and approaches zero from below for 
points far away from the substrate (Fig.lT^. 



FIG. 17. Dependence of the radius R of curvature (a) and of 
the maximal height h (b) of the liquid-gas interface shown in 
Fig. union the reduced volume V — V/{Lyb^) of the droplet 
for 6 > tp. The length b is half of the width of the basis of the 
trapezoid. Vi and V2 are the limiting values corresponding 
to configuration II for which the three-phase contact line is 
pinned to the edge of the substrate so that the apparent con- 
tact angle a varies in the range 9 ^ a ^ 9 -'r ip (see Fig. ll6|) . 
depending on the volume V. Vi* = V{(x = 9)/{Iyb'^) and 
V2* = Via ^ 9 + ip)/{Lyb^) (see Eq. In (6), h/b in- 

creases for / and // whereas the curve for /// increases for 
9 > tp, IS constant for 9 = p>, and decreases for 9 < p>. In 
both (a) and (b) the angles 9 and ip are chosen as 9 = 7r/4 
and ip = 7r/6. 



2. Mesoscopic description 

In the mesoscopic description one takes into account 
the presence of the wetting film the droplet is connected 



3. Equilibrium shape of the sessile droplet 

The effective Hamiltonian for the interface of a sessile 
droplet deposited on a trapezoidal substrate is given by 



Mrap[f]^Ly ( dx[-l[^l + {.f'{x)Y-^l + {{a'{x)y 



+ Wtrap{x, f{x))-UJtrap{.X,a{x))^ . 

(90) 

As the vertical distance from the planar base of the trape- 
zoid the function / — f{x) describes the shape of the 
liquid-gas interface and 

a{x) = i^e(^{b + tan ^)-\x\^ 



(6— \x\) tan ip+- 



COSif 



e(|a;|-(6+£ootan|) 
(91) 
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FIG. 18. Disjoining pressure for a trapezoidal substrate 
with a basis of width 2b = 20ct and a characteristic an- 
gle (p = 7r/6. Positive values of the disjoining pressure are 
summarily marked with red color. The fluid- fluid and the 
substrate-fluid pair potentials are of the van der Waals type 
with parameters rendering a thickness lo — 2a of the wetting 
film on a planar substrate and a contact angle 9 — 7r/4. 



describes a continuous reference configuration (see the 
violet line in Fig.ll9p which contains an a priori unknown 
thickness £oo as a parameter. The excess volume 



dx 



f{x) - a{x) 



(92) 



is fixed. The equilibrium profile / 
equation (see Eq. ([87]) ) 



/(x) satisfies the 



du!trap{x, Z) 



(l + /'(a-)2)3/2 



^tTap{xJ{x)) - 



(93) 



A. 



where A is a Lagrange multiplier. In the following we 
omit the overbar indicating the equilibrium profile. 

The equilibrium shape of the interface is taken to be 
symmetric with respect to the plane a; = together with 
f'{x = 0) = 0. Moreover we assume that for x — ?> ±cxj the 
shape f{x) approaches the function a{x). For x — )■ ±oo 
the effective interface potential converges to its planar 
substrate form 

liin UJtrap{x,f{x)) = UJ.p( COS'f{f{x)-a{x))+£oo) , 

(94) 

so that the Lagrange multiplier A can be deter- 
mined from the boundary conditions at infinity. Since 
f"{x ^ ±oo) = and Up{e) = -w^(^), Eq. leads to 



A = t^:(€oo). 



(95) 



For a given value of £oo the solution f{x) of the 
second order differential equation (|93p contains no free 



parameter. The two integration constants are deter- 
mined by the boundary conditions /'(x = 0) =0 and 
f{x — > — oo) = {x + b) tan ip + laoj cos (/s. On the other 
hand the parameter irx, is determined by the excess vol- 
ume V^r^ of the droplet (Eq. (|92|) ). 

In order to find the equilibrium profile of the liquid-gas 
interface we use a procedure analogous to the one used 
in Sec. lIIIBll 

1. fix /(x = 0) at a certain value /o > 0; 

2. fix irx, > io and calculate the Lagrange multiplier 
A (Eq.®); 



3. integrate Eq. ((93)) numerically with the boundary 
conditions /(O) = /o and /'(O) = within the range 
X £ [—L, 0], where x = —L is the imposed limit of 
the system size on the left hand side; 

4. compare ,f{—L) with a{—L) and f'{—L) with 
a'{—L) = tan ip; 

5. if the differences |/(— L)/a(— L) — 1| and 
\f'{—L)/a'{—L) — 1| arc not satisfyingly small 
return to step 2 and use a different choice for £oo- 

This procedure fixes f{x — 0), and £ao and thus Vex fol- 
low; this relationship can be inverted. Our method is re- 
stricted to contact angles 9 within the range 7r/2 — (p]. 
At the upper limit tt/2 — if the liquid-gas interface can 
develop overhangs which cannot be described by a single- 
valued function / = f{x). On the other hand, for 9 ^ ip 
there arc many loo corresponding to the same /□ and 
it is not obvious how to choose the new value of £ca 
when returning from step 5 to step 2 in the above al- 
gorithm. This problem becomes already apparent within 
the macroscopic description according to which in the 
case 9 < ip the height h of the droplet is the same for two 
distinct volumes (see /// in Fig.[T7](b)). 

For each solution of the equilibrium shape of 
the liquid-gas interface we fit the arc of a circle 



R+ fo with the same curvature as 



the one of the liquid-gas interface at a; = fFig. fTO]) . The 
radius of the arc of the circle is given by (see Eqs. ([93)) 
and dMI) 



R 



7 



n(rap(0, fo) 



(96) 



In addition, for each equilibrium profile the position 
X — Xd < a,t which the arc of the circle intersects 
the asymptote function, ar{xd) = a{xd), and the corre- 
sponding angle a = arctana^(a;d) are determined. Upon 
increasing /o (i.e., increasing the volume of the droplet) 
the position Xd moves smoothly across the edge of the 
trapezoidal substrate and the angle a increases (Fig. [20| . 
The latter is bounded from above hy 9 + p> in accordance 
with Gibbs' criterion. However, for Xd sufficiently far to 
the right of the edge one finds a < 9 (Fig. [^D|) which 
is in contradiction with Gibbs' criterion for macroscopic 
droplets. 
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FIG. 19. Equilibrium shape f{x) of the liquid-gas interface 
for fo = 17, which renders loo = 2.04, and the fitted arc of 
a circle ar{x) for a sessile droplet deposited on a trapezoidal 
substrate with a characteristic angle = 7r/6 and a basis 
width 26 = 100(7. The thickness of the wetting film far away 
from the edge is denoted by £oa ■ The contact angle a is defined 
as the slope of ar{x) at the intersection x = Xd of ar{x) with 
the asymptote a{x) (Eq. H91I) ). The parameters of the effective 
interface potential are aa — 0.5, 6 = n/A, and £o = 2a. 



For the equilibrium solutions characterized by fo (or, 
equivalently, ioo or Vex) we have calculated the contact 
angle 6d{£oo) of the droplet deposited on a planar sub- 
strate (Eq. ([74]) ). It turns out that the function a{xd) 
tends to Odiioo) for Xd + b > (Fig.[20l). This behav- 
ior can be understood by noting that for nanodroplets 
the contact angle 6d changes significantly vi^ith their vol- 
ume (Fig. [15]). The application of Gibbs' criterion also to 
nanodroplets states that the contact angle of the sessile 
droplet near the edge is bounded from below by the con- 
tact angle of the corresponding droplet deposited on the 
same but planar substrate and from above by the contact 
angle 9 + ip, as for macroscopic droplets. 

We recall that the contact angle of nanodroplets de- 
pends sensitively on their volume. For macroscopic 
droplets 9d ^ and in Fig. [20] the difference between 
Young's angle 9 and 9d{ioo) would vanish. In this limit 
the shape of the function a{xd + b) would resemble the 
one obtained for the liquid-gas interface adsorbed at an 
apex-shaped substrate (Fig. [7]). Macroscopic droplets de- 
posited on a trapezoidal substrate with Xd + b > can 
be prepared for macroscopic values of the width 2b. 

For simplicity instead of the excess volume Vex we use 
the volume Vd defined as the excess volume correspond- 
ing to the function ar{x) over the function a{x) in order 
to characterize the volume of the droplet approximately. 
Within a mesoscopic description the radius of curvature 
R, the contact angle a, and the parameter Xd are smooth 
functions of Vd, which is in contrast to the macroscopic 
description (Figs. [2T]-[23]). For large volumes these quan- 
tities approach their macroscopic analogues. It turns out 
that for small volumes they can be described by the cor- 
responding macroscopic equations if one uses the volume 
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FIG. 20. Dependence of the contact angle a on the shifted po- 
sition Xd + b of the three-phase contact line for the equilibrium 
shape of the sessile nanodroplets as the one shown in Fig. 1191 
The horizontal lines denote the angles 6 and 9 + ip expected 
from Gibbs' criterion. The short line oi = Od{£oo) on the right 
end corresponds to the contact angles for nanodroplets de- 
posited on the same but planar substrate for the same values 
of loo as the ones rendering the data points there, for which 
there is a one-to-one correspondence between Xd and loo as 
well as between loo and Vex (see for comparison the lower 
curve in Fig. ll4l corresponding to a planar substrate). The 
parameters of the effective interface potential and the angle 
are the same as in Fig. 1191 and b = 50cr. 



dependent contact angle 9d. 

For effective interface potentials rendering 9 > ip the 
height fo of the droplet is an increasing function of 
its volume Vd fFig.lMTa')') as in the macroscopic case 
(Fig.^nib)). The thickness (Fig-Mib)) is a non- 
monotonic function of the volume (in contrast to the 
planar case (see Fig.[H|). which signals the transition 
from configurations / to /// introduced for macroscopic 
droplets (Scc. lIVBip . 
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FIG. 21. Dependence of the radius of curvature R of liq- 
uid droplets deposited on a trapezoidal substrate (Fig. ll9|) on 
their volume Vd within the mesoscopic (dots) and the macro- 
scopic (full lines, Fig. ll7f a)) description. The line Id denotes 
the radius calculated within the macroscopic description but 
taking into account the change with volume of the contact 
angle 9d of the corresponding nanodroplet on the planar sub- 
strate. The surface tension coefficient, the parameters of the 
effective interface potential, and the angle are the same as 
in Fig. 1191 The droplet volume is measured in units of Lya^. 




FIG. 23. Dependence of the distance Xd + b from the edge of 
the substrate of the three-phase contact line (Fig. ll9|) of liquid 
droplets deposited on a trapezoidal substrate on the droplet 
volume Vd within the mesoscopic (dots) and the macroscopic 
(full lines) description. The full line // describes the macro- 
scopic pinning for Vi <V <V2 (compare Fig. ll7p . The short 
line Id denotes Xd + b as calculated within the macroscopic 
description but taking into account the change of the contact 
angle of the nanodroplet 9d with volume. The parameters of 
the underlying effective interface potential and the angle ip 
are the same as in Fig. 1191 The droplet volume is measured 
in units of LyO^ . 
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FIG. 22. Dependence of the contact angle a of the liquid 
droplet deposited on a trapezoidal substrate on the droplet 
volume Vd within the mesoscopic (dots) and the macroscopic 
(full line) description (compare Fig. [T7)l . The short line dd{(.oo) 
denotes the contact angle calculated for the droplets deposited 
on the same but planar substrate with the values of loo cor- 
responding to the ones for the dotted line. (According to 
Fig. [20] there is a one-to-one correspondence between i^o and 
Vd — Vex-) The parameters of the underlying effective inter- 
face potential and the angle tp axe the same as in Fig. 1191 Fig- 
ure [22] translates Fig.[2S] into the dependence on the droplet 
volume, which is measured in units of LyO^ . 




FIG. 24. Dependence of the maximal height fa (a) and of the 
wetting film thickness i^o (6) of liquid droplets deposited on 
a trapezoidal substrate on their volume Vj. The parameters 
of the underlying effective interface potential and the angle ip 
are the same as in Fig. ll9l 
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4. Width of the transition region 

Within the macroscopic description the three-phase 
contact Une remains pinned at the edge of the trapezoidal 
substrate for a certain range Vi < V < V2 of droplet 
volumes (see Fig. [23| . For such configurations (denoted 
as // in Sec lIVBT]) the radius of the droplet decreases 
with its volume (Fig.[T7Ja)) while for configurations / 
and III it is an increasing function of the volume. In 
the mesoscopic description there is no contact line pin- 
ning. However, the radius R of the arc of the circle fit- 
ted to the equilibrium liquid-gas interface has a similar 
non-monotonic volume dependence as in the macroscopic 
description (Fig.[2l])- According to Figs.[2T] and [23l the 
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FIG. 25. Dependence of the radius R of the arc of the circle 
fitted to the interface shape on the distance Xd + b of the 
three-pliase contact line from the edge of tlie substrate (see 
Fig.ll9[). The distance w between the abscissa of the local 
extrema is a measure of the width of the transition region. 
The parameters of the underlying effective interface potential 
and the angle are the same as in Fig. 1191 



region where R is an increasing function of -I- ^ corre- 
sponds to the aforementioned pinning in the macroscopic 
description. The spatial extent of this region is denoted 
by w and is a measure of the width of the transition re- 
gion within which the contact line passes smoothly across 
the edge (Fig. [25]). For the choice of the parameters used 
in Fig. [^S] the width of the transition region is of the order 
of ten fluid particle diameters and thus is mesoscopic in 
character. 



5. Line contribution to the free energy 

For the system under investigation we define the line 
contribution to the free energy as the difference between 
the free energy of the droplet and the free energy of the 
configuration corresponding to the arc of the circle ar{x) 
fitted to the droplet, both relative to the reference con- 



figuration a(x): 

T[f] —{■y^trapif] — -J^traplar]) / Ly = 

2 r dx + (/'(x))2 - Vl + (a'(x))2" 

J — 00 

(x, f{x)) (a;,a(a;))| 



2 / da 



\x {7 [Vf + «(x))2 - + {a'{x)f 

+ UJtrap{x,ar{x)) - LOtTap{x , a{x))^ . 



7(7ct) 



(97) 
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FIG. 26. Dependence of the line contribution r on the dis- 
tance a;o -f 6 of the three-phase contact line from the edge of 
the substrate of a sessile droplet deposited on a trapezoidal 
substrate. The parameters of the underlying effective inter- 
face potential and the angle tp axe the same as on Fig. 1191 



The line contribution changes significantly when the 
three-phase contact line passes the edge (Fig. [26]). The 
spatial extent of the region in which the line contribution 
is an increasing function of the volume of the droplet (i.e., 
decreasing as function of Xd + h) is of the order of three 
fiuid particle diameters. If the contact line is far from 
the edge the line contribution is a decreasing function of 
the volume of the droplet. 

We emphasize that the line contribution to the free 
energy presented in Fig. [25] corresponds to different equi- 
librium profiles of the droplets, in particular with differ- 
ent volumes. It is not the plot of the line energy of a 
droplet with a fixed excess volume. Nonetheless, Fig. [26] 
indicates that there is a free energy barrier at the edge of 
the substrate and thus a moving droplet with fixed excess 
volume is expected to stop just before reaching i t^^'^° . 



V. SUMIVIARY AND DISCUSSION 
A. Summary 

If a substrate surface forms a sharp corner and 
the three-phase solid-liquid-gas contact line of a sessile 
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droplet is pinned at the substrate apex the modified 
Young's equation for the contact angle (Eq. ([!])) is no 
longer valid. Instead, the corresponding local contact an- 
gle a can take any value within the range 9 ^ a ^ 6 + tp 
(Eq. (m), where n — ip is the opening angle of the apex 
formed by the substrate faces (Fig.[T]). This ambiguity 
of the local contact angle at the apex is called Gibbs' 
criterion. In order to determine the equilibrium shape 
of the liquid-gas interface of a liquid film covering an 
apex-shaped substrate (Fig. [5]) and the equilibrium shape 
of a sessile droplet deposited on a trapezoidal substrate 
fFig.irS]) we have used an effective interface Hamiltonian 
based on density functional theory. This approach has 
proved to be very useful in analyzing similar systems^. 
The thermodynamic state of the system is taken to be 
at the bulk liquid-gas coexistence line below the wetting 
temperature and well below the critical point of the liq- 
uid. For our explicit calculations we have chosen the 
thickness of the wetting film to be of the order of a 
few fluid particle diameters. We have focused on quasi 
two-dimensional systems being translationally invariant 
in one direction. 

First, we have analyzed the equilibrium shape 
f{x) of the liquid-gas interface at a planar substrate 
(Fig. 131) with boundary conditions f{x — > — oo) = io and 
/'(x — >■ (X)) = tan^. The film thickness £o minimizes 
the effective interface potential for a planar substrate 
uip{£) and the angle 6 fulfills the macroscopic Young's 
law (Eqs. (HH) and ([20])). If the height of the interface 
at one point xq is fixed as /o = /(xq), the derivative 
/o ~ f'i^o) is a unique function of /o (Fig.S]) for the afor- 
mentioned boundary conditions f{x — >■ — oo) = io and 
/'(a; — > oo) = tan(?. It changes monotonously from 
/o = for /o = £o to /o = tanO for /o oo. 

In the case of an apex-shaped substrate we have calcu- 
lated the equilibrium profile for the liquid-gas interface 
numerically. Due to limited numerical accuracy, for a 
fixed value /g one cannot find the value of the derivative 
/o rendering the exact boundary condition on the far left 
hand side x = Li of the system. What can be achieved 
numerically is to find values /q = and /q = ren- 
dering solutions f<{x) and f>{x), respectively, which for 
sufficiently large, negative x follow the asymptote aap{x) 
(Eq. (j46p ) and differ only slightly from it in the vicinity 
of the boundary x = Li (Fig. [5]). The equilibrium pro- 
file f{x) lies between these functions f<{x) and fy{x). 
For each solution f{x) we determine the contact angle 
a = arctan/'(L2) (where L2 is the imposed boundary at 
the right hand side of the system) and the quantity Xd 
characterizing the position at which the liquid-gas inter- 
face detaches from the substrate (Fig. [6]). The contact 
angle a is a decreasing and continuous function of Xd 
(Fig. [7]); there is no indication of three-phase contact line 
pinning. For Xd — ±00 the contact angle a approaches 
from below its limiting values 9 =F </3/2, which arc those 
expected from Gibbs' criterion for this geometry. 

In Sec. IV we have examined cylindrical droplets with 
a fixed excess volume as well as incomplete wetting films 



deposited on planar and trapezoidal substrates. In the 
planar case (Fig. [5]) we have obtained the equation for 
the shape of the interface for various excess volumes of 
the droplel42ri&i2£. Both in a finite system of width 2L 
and in the case when the droplets can spread over an un- 
bounded substrate, i.e., L — >■ 00, characteristic features 
of the droplet shape can be inferred from tangential con- 
structions to the effective interface potential u)p{t) (Figs. 
l9landfT2)). This involves in particular the heights /l (i.e., 
H.00 in unbounded systems) and /o denoting the minimal 
and maximal values, respectively, of the function f(x) 
describing the shape of the equilibrium nanodroplets. In 
laterally unbounded systems, /o — >■ 00 and — >■ ^0 for 
increasing excess volumes Vex ~^ 00 fFig.lH]). 

For a laterally finite, planar system of size L the 
droplet solution does not exist for arbitrary > £q. 
If /l is sufficiently small there is no /o which would cor- 
respond to the prescribed fixed L (Fig. [TU|) . On the other 
hand the quantity Jl is bounded from above, rendering 
upper (Vt™''^(L)) and lower (\/™"(i)) bounds for the 
total volume for which the droplet configuration exists 
(Fig.HH). 

We have also calculated the contact angle for nan- 
odroplets (Fig. [131). It is smaller than Young's angle for 
macroscopic dropletSii^ and it is an increasing function 
of the volume of the droplet (Fig. [15]). This dependence 
has to be taken into account also for analyzing sessile 
nanodroplets deposited on trapezoidal substrates. 

For symmetrical ridgelike macroscopic droplets de- 
posited on trapezoidal substrates one can distinguish 
three different configurations depending on the position 
of three phase contact line (Fig. fTB)) . The radius R of 
curvature of the droplets is a continuous but neither 
a smooth nor a monotonous function of the volume of 
the droplet. It is decreasing if the three-phase conctact 
line is pinned to the edge of the substrate (configuration 
//), and increasing otherwise (configurations / and ///) 
(Fig.[I7D. 

In order to determine the equlibrium shape of nan- 
odroplets on trapezoidal substrates we have calculated 
the disjoining pressure for such substrates fFig. fTS)) . For 
each ensuing equilibrium shape of the liquid-gas inter- 
face we have fitted the arc of a circle with the same ra- 
dius of curvature R as the one of the liquid-gas interface 
at the center x = (Fig.^^. In addition, the position 
X = Xd < &t which the arc of the circle ar{x) inter- 
sects the asymptote a{x) of the film thickness and the 
corresponding contact angle a = arctan a[.(xd) were de- 
termined. Upon increasing /o (i.e., increasing the volume 
of the droplet) the position Xd, which can be interpreted 
as the three-phase contact line position, moves smoothly 
across the edge of the trapezoidal substrate and the con- 
tact angle a increases fFig. l^O)) . The latter is bounded 
from above by 9 + ip in accordance with Gibbs' criterion. 
However, for Xd sufficiently far to the right of the edge 
one finds a < 9 which is in contradiction with Gibbs' 
criterion for macroscopic droplets. This behavior can be 
understood by noting that for nanodroplets the contact 
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angle 6d changes significantly with their volume (Fig.[T5|). 
The extension of Gibbs' criterion to nanodroplets states 
that the contact angle of sessile droplets near the edge 
of the substrate is bounded from below by the contact 
angle of the corresponding finite-sized droplets deposited 
on the same but planar substrate and from above, as for 
macroscopic droplets, by the contact angle 9 + ip. 

Within a mesoscopic description the radius of curva- 
ture R, the contact angle a, and the parameter Xd arc 
smooth functions of the volume of the droplet Vd, which 
is in contrast to the macroscopic description (Figs.[5T]- 
[^5)1 . For large volumes these quantities approach their 
macroscopic analogues. It turns out that in the limit of 
small volumes they can be described by the correspond- 
ing macroscopic equations if one uses the volume depen- 
dent contact angle 6d of nanodroplets. 

For effective interface potentials rendering 9 > (p the 
height /o of the droplet is an increasing function of its vol- 
ume Vd (Fig.[24|). as in the macroscopic case (Fig. [17^6)). 
Opposite to the planar case fFig.lH)) the film thick- 
ness £oo is a non-monotonic function of the drop vol- 
ume, which signals the transition from configuration / to 
/// introduced for macroscopic droplets. According to 
Fig. [21] the region where i? is a decreasing function of Vd 
corresponds to three-phase contact line pinning within 
the macroscopic description (configuration //). The spa- 
tial extent w of this region is a measure of the width of 
the transition region within which the contact line passes 
smoothly across the edge (Fig. [25]) . For the choice of the 
parameters used in Fig. [55] the width of the transition re- 
gion is of the order of ten fluid particle diameters and 
thus is mesoscopic in character. 

The line contribution to the free energy of the droplet 
changes significantly when the three-phase contact line 
passes the edge (Fig. [55]). The spatial extent of the region 
in which the line contribution is an increasing function of 
the volume of the droplet (i.e., decreasing as function of 
the distance Xd + b oi the three-phase contact line from 
the edge) is of the order of three fluid particle diameters 
and thus also mesoscopic in character. If the contact line 
is far from the edge the line contribution is a decreasing 
function of the volume of the droplet. The edge poses a 
free-energy barrier for the three-phase contact line. 

Our numerical results for Gibbs' criterion have been 
obtained from the analysis of the disjoining pressure for 
apex-shaped substrates for specific choices of the effective 
interface potential (based on long-ranged interparticle in- 
teractions) and for specific geometrical parameters. We 
have studied the absence of three-phase contact line pin- 
ning on the nanoscale and we have analyzed how Gibbs' 
criterion has to be modified in order to describe sessile 
nanodroplets on substrates with sharp edges. 



B. Discussion 

Three-phase contact line pinning which takes place at 
asperities of non-planar, chemically homogeneous sur- 



faces is a common phenomenon due to inherent rough- 
ness of both naturally occuring as well as fabricated sub- 
strates. Important examples vary from capillary filling of 
geometrically patterned channels to terraced substrates. 
In the case of channels patterned by pillars, depending 
on the distance between these obstacles, the width of the 
channel, and Young's contact angle the advancing liquid 
front can be pinned and flow can be suppressedi^^i^. In 
Sec. lIIII we have shown that within a mesoscopic descrip- 
tion there is no three-phase contact line pinning of the 
liquid-gas interface at a substrate edge due to the extra 
cost related to the associated increase of the line contri- 
bution to free energy of the system. Thus one can specu- 
late that dense arrangements of obstacles favor capillary 
flow in microchannels despite the limitations predicted 
by the macroscopic version of Gibbs' criterion. This con- 
jecture seems to be supported by the fact that dewctting 
of terraced substrates can proceed for step heights up to 
a couple of nanometers while under the same thermody- 
namic conditions dewetting is suppressed for larger step 
heights^. 

Three-phase contact line pinning has also an impact 
on contact angle hysteresis on rough surfacesi^. Be- 
sides the shape and the chemical character of the as- 
perities, their size plays an important role. Already 
in early studies of three-phase contact line pinning, it 
was shown that steps with a height of approximately 
10 nm do not pin the three-phase contact line at the 
edge of stepped substrate^. Advances in the fabrica- 
tion of nanostructured substrates have allowed the inves- 
tigation of the contact line behavior at rings grown on 
a flat substrate with a trapezoidal vertical cross section. 
For rings with heights below 2 the advancing contact 
angle decreases signiflcantly with the height of the trape- 
zoidal asperity2^. Both these observations deviate from 
the macroscopic Gibbs' criterion and can be related to 
the findings presented in Sec. lIVi according to which the 
position of the three-phase contact line moves continu- 
ously with a smooth variation of the apparent contact 
angle near the edge of the substrate. More quantitative 
analyses concerning the suppression of three-phase con- 
tact line pinning at nanometer sized steps are warranted 
and promising. 

For sessile nanodroplets with a fixed excess volume de- 
posited on a laterally unbounded apex-shaped substrate, 
there arc three morphologically distinct solutions of the 
Euler-Lagrange equation (Fig. l27|) . Two of them are sym- 
metric. For the same excess volume the three configura- 
tions have different asymptotic thicknesses £oo and there- 
fore it is not clear a priori which configuration has the 
lowest free energy and thus is the stable one. 

The disjoining pressure for the liquid-gas interface at 
a trapezoidal substrate can be expressed in terms of the 
difference of the disjoining pressures stemming from two 
suitable apex-shaped substrates (Eq. ([87]) '). Thus also for 
a trapezoidal substrate there are different configurations 
of the liquid-gas interface fulfilling the Eulcr-Lagrange 
equation, as for the apex-shaped substrate case. 
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(a) 




FIG. 27. Schematic shapes of three morphologically distinct 
equilibrium droplet configurations with the same fixed excess 
volume deposited on an apex-shaped substrate. The configu- 
rations (a) and (b) are symmetric with respect to the mid- 
plane. The three configurations have different asymptotic 
thicknesses: ^1^^ , l'^ , and . 



Here our investigation has been focused on sessile 
nanodroplets which are symmetric and attain their 
maximal height at the center of the system. The issue 
of morphological transitions between different sessile 
droplet configurations on apex-shaped and trapezoidal 
substrates is left for future research. This has been 
already investigated for macroscopic droplets deposited 
on axisymmetric pillar substrates^ii^. Recent studies of 
the Vapor-Liquid-Solid mechanism of nanowire growth 
show that the liquid droplet promoting the solid growth 
can wet the sidewall of the nanowire and thus does not 
sit at the top of the pillar with the three-phase contact 
line pinned to its edge as for typical VLS growth^. 
A theoretical description of the transition between 
these two configurations using the present mesoscopic 
approach appears to be interesting. 

Finally, we mention an interesting process in which 
a droplet of fixed volume V is placed on a trapezoidal 
substrate and its contact angle 9 is decreased, e.g., by in- 
creasing an applied voltage as in electrowetting22r2^. Ini- 
tially, the droplet shape corresponds to configuration / 
on Fig.[Tni Upon decreasing the contact angle the droplet 
spreads until the three-phase contact line reaches the 
edge of the trapezoidal substrate. If the corresponding 



contact angle fulfills 9 = a > ip (where a is the solution of 
Eq. (HI])), upon further increase of the voltage, the shape 
of the droplet and the apparent contact angle remain con- 
stant until the angle decreases to the value a — ip (in 
agreement with Gibbs' criterion). Upon further increase 
of the voltage one expects the contact angle 9 to start 
to decrease again and the drop to spread on the tilted 
side of the trapezoidal substrate. In actual experimental 
settings the above naive scenario may be substantially 
modified by effects related to the fact that the drop is 
charged. In particular the change of morphology of the 
droplet front upon reaching the apex in the presence of 
electric fields provides interesting scientific perspectives. 

In summary, we have shown that the presence of meso- 
scopic wetting films on edged substrate surfaces prevents 
three-phase contact line pinning on the nanoscale. We 
have analyzed the shape of the liquid-gas interface of liq- 
uid films both at an apex-shaped substrate and for liquid 
sessile nanodroplets with fixed excess volume deposited 
on trapezoidal substrates. Near the edge of an apex- 
shaped substrate the apparent contact angle changes con- 
tinuously within the range of values expected from Gibbs' 
criterion while the three-phase contact line smoothly 
passes through the atomically sharp apex. For a trape- 
zoidal substrate, upon increasing the volume of the nan- 
odroplet the apparent contact angle fulfills a modified 
Gibbs' criterion, for which one has to take into account 
the dependence of the contact angle of the nanodroplet 
on its volume. In both cases the spatial extent of the 
region, within which the three-phase contact line passes 
across the edge, is of the order of ten fluid particle diam- 
eters and thus is mesoscopic in character. 
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